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ABSTRACT 

AE Aquarii is a propeller system. It has the shortest spin period among cataclysmic 
variables, and this is increasing on a 10^ yr timescale. Its UV spectrum shows very 
strong carbon depletion vs nitrogen and its secondary mass indicates a star far from 
the zero-age main sequence. We show that these properties strongly suggest that 
AE Aqr has descended from a supersoft X-ray binary. We calculate the evolution 
of systems descending through this channel, and show that many of them end as 
AM CVn systems. The short spindown timescale of AE Aqr requires a high birthrate 
for such systems, implying that a substantial fraction of cataclysmic variables must 
have formed in this way. A simple estimate suggests that this fraction could be of order 
one-third of current CVs. We emphasize the importance of measurements of the C/N 
abundance ratio in CVs, particularly via the Civ 1550/Nv 1238 ratio, in determining 
how large the observed fraction is. 

Key words: novae, cataclysmic variables — binaries: close — stars: evolution — 
stars: individual: AE Aqr — abundances — stars: individual: AM CVn 



1 INTRODUCTION 

AE Aqr is one of the most distinctive cataclysmic variables 
(CVs). It has a long orbital period of 9.88 hr (Welsh et al., 
1993), and the shortest coherent pulse period (33 s), increas- 
ing on a timescale ^ 10^ yr (de Jager et al., 1994). Doppler 
tomography reveals the apparent absence of an accretion 
disc. This has led to its interpretation as a 'propeller' system 
(Wynn et al., 1997): the white dwarf is apparently expelling 
centrifugally the matter transferred from the companion. To 
arrive at such a highly non-equilibrium spin rate, the mass 
transfer rate in the recent (lO'^ yr) past must have been much 
higher than its current value M2 ^ few x 10~^ MQyr~^, and 
decreased on a still shorter timescale. lUE spectra of AE Aqr 
(Jameson et al., 1980) reveal the most extreme Civ to Nv 
ratio of all CVs (Mauche et al., 1997), probably indicating 
strong carbon depletion and thus CNO cycling. The short 
spin period makes the system effectively a double-lined spec- 
troscopic binary, with masses M\ ~ 0.89 ± 0.23 Mq for the 
white dwarf (WD) and M2 ~ 0.57 ± 0.15 for the donor 
star (Welsh et al., 1995; Casares et al., 1996). The donor 
star is of spectral type K5V (Welsh, 1999) which is too late 
for the given orbital period and too early for the measured 
M2 if it is on the main sequence. 

We show here that all of these properties are consistent 



with the idea that AE Aqr descends from a supersoft X-ray 
binary, cf. Schenker & King (2002). The thermal-timescale 
mass transfer characterizing such sources ends once the sec- 
ondary/primary mass ratio q = M2/M1 decreases suffi- 
ciently, leading to a rapid transition to normal CV evolution 
driven by angular momentum losses, cf. King (1988). If the 
white dwarf is magnetic, as in AE Aqr, the high mass trans- 
fer rate in the thermal-timescale phase will have spun it up 
to a short spin period. This must lengthen rapidly in order 
to reach equilibrium with the lower transfer rate in the CV 
state, and this spindown is what drives its propeller action. 
The very short lifetime of the propeller phase implies a high 
birthrate for such systems, comparable to those of known 
CVs. This in turn suggests that systems descending from 
evolutions like AE Aqr must constitute a large fraction of 
all CVs. 

The simple idea that all CVs form with essentially un- 
evolved, low-mass main sequence (MS) donors has been 
repeatedly challenged for more than a decade (Pylyser & 
Savonije, 1988a,b; Baraffe & Kolb, 2000; Schenker, 2001; 
King & Schenker, 2002; Schenker & King, 2002). If mass 
transfer starts with mass ratio g ^ 1 the secondary's Roche 
lobe -Rl tends to shrink with respect to its thermal- 
equilibrium radius -Rte- Mass is therefore transferred on a 
thermal timescale (~ few xlO^ yr), at rates high enough 
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to sustain steady nuclear burning on the material accreted 
by the white dwarf, and thus explain (short-period) super- 
soft X-ray sources (van den Heuvel et al., 1992; King et al., 
2001). 

After this relatively brief phase, q becomes small enough 
that i?L shrinks more slowly than Rte, and thermal- 
timescale mass transfer ends. The systems either evolve off 
to longer orbital periods, driven by nuclear evolution, or 
switch to stable, angular-momentum-loss-driven mass trans- 
fer like ordinary CVs. We will show that AE Aqr is at the end 
of this transition phase to a CV, accounting for its strange 
properties. AE Aqr is passing through this transition very 
rapidly; its whole evolution up to its current state is short 
compared to the stable CV phase afterwards. In this sense 
AE Aqr is not unique at all: a significant fraction of CVs 
must descend from similar evolutions. 

In the next section we derive some rough estimates for 
the potential progenitor system of AE Aqr, using a sim- 
plified analytic description of the orbital evolution during 
rapid mass transfer. Sect. 3 studies the behaviour of sin- 
gle star models under fixed mass loss of the order expected 
during such an evolution. This allows us to identify the in- 
fluence of varying stellar parameters, and establish whether 
the observed properties of AE Aqr are compatible with this 
scenario. We present a set of full binary evolution calcula- 
tions in Sect. 4 and discuss them in Section 5. Section 6 is 
the Conclusion. 



2 SIMPLIFIED ORBITAL EVOLUTION 

To get a first idea about the required initial parameters 
for the progenitor system of AE Aqr, we recall that during 
thermal-timescale mass transfer additional angular momen- 
tum losses (magnetic braking & gravitational radiation) are 
negligible compared to the changes of the Roche radius re- 
sulting from mass transfer itself, e.g. King et al. (2001). We 
define tj by 



Ml = -rjM2 



(1) 



as the fraction of the mass lost from the secondary that is 
actually accreted on to the WD, while assuming that the 
rest leaves the system with average specific angular momen- 
tum equal to the primary's, i.e. /3 = 1 in the nomenclature 
of King et al. (2001). From their eq. (10) we obtain the evo- 
lutionary tracks as 
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respectively for > 0. For the limit = these lead to the 
exponential expressions given e.g. in King & Ritter (1999). 

Equipped with these expressions we can estimate where 
analytically computed curves for the evolution of radius and 
orbital period intersect those describing mass loss under the 
assumption of permanent thermal equilibrium (i.e. essen- 
tially along the MS). These intersections mark the transi- 
tion between the fast thermal-timescale mass transfer (su- 
persoft) and CV-like angular-momentum-loss-driven evo- 
lution at lower mass transfer rates. Figure 1 shows a set 




Figure 1. Evolution of a model AE Aqr progenitor system. For 
the assumed parameters and A/i^now = 0.89 Mq, various tracks 
are shown, starting from a 1.6 Mq star which has almost reached 
its maximal MS radius. The three panels show the evolution of the 
orbital period, secondary radius and WD mass for different cases 
(dotted line; r) = 1 — full line: rj = — dashed line: rj = 0.3). 
The diamond marks the current position of AE Aqr. Note that 
the conservative model {rj = 1) cannot meet the requirements 
for Mi: even starting from the lowest plausible WD mass, it has 
grown well beyond the Chandrasekhar limit before reaching the 
current Porb of AE Aqr. The grey shaded area in the middle panel 
marks the radii of single stars during their main sequence life as 
labelled. 



of curves for different 77 together with the shaded location 
of the MS framed by a pair of thick curves marking the 
zero-age MS (ZAMS) and the turn-off MS (TMS). For the 
latter we sue the point of maximal radius prior to core hy- 
drogen exhaustion. Note that this definition is not applica- 
ble to MS stars without convective cores (i.e. with masses 
below ~ 1.2 Mq). The lowest panel shows the growth of 
the WD mass as the donor evolves (from right to left) 
for various rj. The initial WD mass Mi, is chosen so that 
the current position of AE Aqr is reached with masses of 
(Ml, Ma) = (0.89Mq,0.57Mq) in accordance with observa- 
tion. A good fit should also pass through the correct orbital 
period and radius (estimated by the corresponding period 
- mean density relation). The assumption of fully conserva- 
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tive mass transfer (77 — 1) leads to much higher WD masses 
than observed, and we therefore do not consider it further. 

Using the curves shown we derive first guesses for the 
system parameters of the progenitor system, i.e. Mu, M2i, 
and ri. We cannot expect these estimates to be very pre- 
cise, as the actual reaction of the system to mass loss in this 
transition phase follows neither the Roche curves nor the 
equilibrium curves (the secondary is out of thermal equilib- 
rium) . 

Nevertheless we can estimate M2i ~ 2 Mq and 
1.5 < Qi < 3, where the upper constraint derives from the oc- 
curence of delayed dynamical instabilities (DDI, see King 
et al. (2001) and discussion section). As explained above, 
conservative mass transfer gives too high a WD mass. How- 
ever a value i] ~ 0.3 is plausible and gives an Mi mildly 
increased over the average field WD mass, in accordance 
with observations of AE Aqr. 



3 SINGLE STAR CALCULATIONS 

We have shown that in principle a recent thermal-timescale 
mass transfer supersoft phase could explain AE Aqr. We now 
use full stellar models. This determines the star's reaction 
to mass loss, and thus whether the observed spectral type 
(SpT) may be reached in this way. We simplify the mass 
transfer history as 

f Ma/TKH VM2 > Mew 

-M2 = < Q 1 (4) 

[2 X lO-^Moyr-i VM2 < M^w 

where tkh is the secondary star's Kelvin-Helmholtz time, 
thus mimicking thermal-timescale mass transfer and the 
subsequent transition to an evolution driven by angular mo- 
mentum loss, as in ordinary CVs. The switching mass A/sw 
is chosen so that thermal-timescale mass transfer ends when 
q — I, i.e. it can be roughly identified with Mi. A more thor- 
ough discussion of the exact stability conditions in CVs is 
presented in Webbink (1985); Sobermann et al. (1997). 

The code we use in both the single star case of the last 
Section and the binary calculations reported here is based on 
Mazzitelh (1989), as adapted by Kolb & Ritter (1992) plus 
minor improvements. The mapping from effective tempera- 
ture to SpT follows the procedure outlined in Beuermann 
et al. (1998) and Baraffe & Kolb (2000). This is not entirely 
self-consistent as our code still uses grey atmospheres. 

At first we performed a number of calculations for Pop. I 
stars with 1.2, 1.6 and 2.0 M0 and initial ages of 0, 0.25, 0.5, 
0.75 and 1 of their MS-lifetime, using a switching mass of 
Msw = 1 Mq. As shown in Fig. 2, more massive stars of the 
same evolutionary phase and further evolved stars of the 
same mass will both lead to larger equilibrium radii after 
the thermal-timescale mass transfer. 

Single star radii can easily be converted into orbital pe- 
riods by equating to the corresponding Roche radius (with 
a primary mass given by Mi). This allows us to com- 
pare our tracks with the observed parameters for AE Aqr, 
Porb = 9.88 hr, SpT = K5 and M2 = 0.57 M©. For the re- 
maining calculations in this section (shown in Figs. 3 and 4) 
a switching mass of 0.7 Mq was used. We varied this value 
slightly for the purpose of testing, but it turned out to have 
little influence on the fitting of final system parameters. 

We obtained four approximate solutions for the initial 




0.0 0.5 1.0 1.5 2.0 

Ms / Me 



Figure 2. Single star calculations with a mass loss rate according 
to eq. (4). The upper panel shows the radius evolution of mass- 
losing stars with initially 1.2, 1.6, 2.0 Mq and an initial age 0.75 
of their MS lifetime. The lower panel illustrates the variation 
with age for a 1.6 Mq star, covering an initial age of 0, 0.25, 
0.5, 0.75 and 1 MS-timescale from bottom to top. Note that the 
strong deviation of the most evolved of the five tracks results 
from ongoing nuclear evolution, which is fast enough because the 
system was already close to entering the Hertzsprung gap. The 
grey shaded area and additional lines mark the location of the 
MS (as labelled). Thermal timescale mass transfer was switched 
off at M2 = 1.0 Mq in each of the computations. 



parameters of AE Aqr (M2,/Mq, Xd, M2/Mq): SI (1.8, 
0.18, 0.6), S2 (1.9, 0.21, 0.6), S3 (2.1, 0.23, 0.54) and S4 (2.3, 
0.247, 0.5), which provided the starting point for the full 
binary calculations in the next section. A decreasing amount 
of initial core hydrogen fraction Xd indicates a donor that is 
advanced futher in its MS evolution. The tracks for S3 and 
S4 are shown in Figs. 3 and 4 as the evolution of spectral 
type over secondary mass and orbital period respectively. 

The (admittedly crude) modelling presented in this sec- 
tion leads to an overall fitting end state of the system, i.e. a 
reasonable model for AE Aqr today. However the preceding 
transition phase is poorly modelled. As a general trend de- 
rived from comparing S1-S4, lower initial secondary masses 
require the donor to be a bit more evolved, and lead to a 
slightly higher current donor mass in AE Aqr. We can there- 



4 K. Schenker et al. 




Figure 3. Evolution of effective temperature (spectral type) with 
decreasing mass in the two best— fitting single star models S3 with 
M2i = 2.1 Mq (dashed line) and S4 with 2.3 (full line). De- 
tails of the conversion from T^ff to SpT are discussed in the text. 
Thermal timescale mass transfer was switched off at M2 = 0.7 Mq 
in accord with the actual WD mass in AE Aqr. Tracks end close 
to the current position of AE Aqr, marked by a diamond. Note 
that the direction of evolution of single stars on the MS (shaded 
area) is towards lower T^ff in more massive stars, but towards 
higher temperatures at lower masses. 
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Figure 4. Evolution of effective temperature (spectral type) with 
orbital period, for the same two tracks S3 & S4 as in the previous 
figure. The large wiggles in orbital period during the early stage 
(at a SpT slightly above GO) are caused by the inconsistency in- 
troduced by an imposed mass transfer rate eq. (4), and the smaller 
ones near SpT KO are directly connected with the jump in M2. 
The exact MS location (shaded area) in orbital period is of course 
a function of the primary mass, for which we have assumed the 
canonical value of Afi = 0.89 Mq in this figure. 



Table 1. Summary of system parameters for selected full binary 
calculation. 



Name 


Mii/M0 


Af2i/M0 




P,/hr 


V 


Bl 


0.95 


2.35 


0.287 


17.9 





B2 


0.95 


2.1 


0.223 


18.5 





B3 


0.95 


1.8 


0.190 


18.1 





B4 


0.8 


2.35 


0.287 


17.5 


0.1 


B5 


0.6 


1.6 


0.560 


10.8 


0.3 


B6 


0.6 


1.6 


0.087 


18.6 


0.3 



fore expect initial values around M2 ~ 2 Mq to provide pro- 
genitors for AE Aqr. 



4 FULL BINARY CALCULATIONS 

There are two major effects suppressed by the single star 
treatment in the previous section: 

(i) For identical initial donor stars, different initial WD 
masses and thus initial mass ratios can produce very differ- 
ent mass transfer rate histories, which are represented too 
simplistically by eq.(4). 



(ii) The growth of the WD mass is ignored in converting 
radii to orbital periods. 

Therefore we have computed a range of tracks in a more 
self-consistent way. The mass transfer rate is calculated ac- 
cording to the current orbital parameters ('full binary evolu- 
tion'). The code used for this work is still based on the Kolb 
& Ritter (1992) version of Mazzitelli's stellar evolution. The 
main challenge in computing the evolutions through and be- 
yond a phase of thermal-timescale mass transfer lies in a 
proper treatment of numerical and physical instabilities. 

We present here the results for six tracks, whose param- 
eters are given in Table 1. In all cases the mass accretion 
efficiency -q defined in eq.(l) is kept constant throughout 
the evolution. Conservative mass transfer (r; = 1) is not 
considered (see Sect. 2). Figure 5 illustrates the variety of 
post-thermal-timescale mass transfer transitions mentioned 
in the discussion of the single star results above. Track B3 
shows just a small bump of enhanced mass transfer. In- 
creasingly S-shaped period loops appear in B2, B6, and Bl. 
Tracks B5 and B4 show the beginning of an apparent run- 
away caused by a delayed dynamical instability (DDI) ; here 
extreme mass ratios coincide with sufficiently deep transient 
outer convection zones. None of these cases, including the 
dip in the mass transfer rate at the end of the thermal- 
timescale mass transfer found in Bl, is adequately described 
by eq.(4). 

The dominant parameter governing the behaviour along 
the sequence B3-B2-B6-B1-B5-B4 is the initial mass ratio: 
larger values correspond to stronger thermal instability. In 
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Figure 5. Possible transitions at tlic end of thcrmal-timcscalc mass transfer shown as the evolution of mass transfer rate over orbital 
period. Each panel is labelled according to system parameters given in Table 1. The features encountered include a simple drop in mass 
transfer rate (B3), a nose-like overhang (e.g. B6), up to cases suffering from an incipient delayed dynamical instability (DDI), countered 
just in time by the reversal of the mass ratio in one case (B5). 



addition the choice of rj (modifying the stability criterion) 
and the evolutionary state of the donor (e.g. B5-B6) can 
also have a strong impact on the evolutionary curves. 

In some of the curves, particularly Bl, a significant 
amount of numerical noise is apparent during the early high 
mass transfer rate phase. This is connected with the growth 
of an outer convcctivc envelope caused by the rapid mass 
loss. As donors of these masses have initially radiative en- 
velopes, very thin convective layers are present at the onset 
of this transition. This creates substantial numerical difficul- 
ties as minute discontinuities in the evolution of the star's 
radius will be reflected in erratic jumps of the mass trans- 
fer rate. As time proceeds, the average thermal reaction of 
the star can be seen to pass right through the average of the 
scatter in A/2: various tests with enhanced resolution and re- 
duced timesteps have confirmed that the subsequent struc- 
ture of the star is unaffected from these problems, which are 
localized in layers close to the surface and thus have quite 
short thermal relaxation timescales. 

Here we discuss tracks B1-B3 and B6. Figure 6 shows 
remarkable similarities to Fig. 3, except for the high ini- 
tial mass ratio cases (Bl & B6) which evolve in a much 
wider arc back to the mass and SpT of AE Aqr. Inset (a) 
shows a zoomed part of the end region of the tracks (dis- 
cussed in the next section), while inset (b) focusses on the 
region near the current position of AE Aqr (again marked 



by an open diamond). Asterisks mark the position along 
each track where the observed mass ratio (presumably the 
best determined observable quantity besides orbital period) 
is matched. Given the uncertainties of spectral type calibra- 
tion in the models and observational mass determinations, 
all four models qualify as potential progenitors. Not supris- 
ingly the evolution vs orbital period shown in Fig. 7 has 
very different features from the corresponding Fig. 4: the 
loops and wiggles in the full binary plot correspond to the 
S-shaped curves shown in Fig. 5. Insets (a) and (b) again 
zoom in on regions of special interest, showing the rather 
good agreement with the observed values of AE Aqr. The 
jitter visible at the beginning of especially B2 results as dis- 
cussed above from numerical difficulties at the onset of mass 
transfer from an outer convective envelope. Again the strong 
influence of a large initial mass ratio on the curves of Bl and 
B6 can be seen. 

All tracks appear very non-MS like after passing 
through AE Aqr, which is a consequence of their advanced 

nuclear evolution. As listed in Table 1 and obvious from 
Fig. 7, the remaining core hydrogen content at the onset of 
mass transfer decreases with initial donor mass (Bl -B2-B3- 
B6) in order to get near the current position of AE Aqr. 

The full binary calculations in this section confirm the 
single star results for a range of suitable progenitors for 
AE Aqr. Evolutionary tracks for a variety of models do pass 
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Figure 6. Tracks of models Bl (dotted), B2 (full), B3 (dash- 
dotted), and B6 (dashed) in the plane of effective temperature 
and secondary mass. The diamond marks the position of AE Aqr, 
the shaded area framed by thick lines indicates the MS. The insets 
provide a zoomed view of the end stages (a) and AE Aqr's position 
today (b). In the latter, asterisks are used to mark the position 
along each track where the observed current mass ratio of AE Aqr 
is reached, showing a fair agreement between the system and the 
proposed progenitor models. 



Figure 7. Tracks of models Bl (dotted), B2 (full), B3 (dash- 
dotted), and B6 (dashed) in the plane of effective temperature 
and orbital period. As in Fig. 4 above, the position of AE Aqr 
and the location of the MS are indicated. Similarly insets provide 
a zoomed view of the end stages (a) and the AE Aqr position 
today (b), in which once again asterisks (marking the observed 
current mass ratio of AE Aqr) confirm the successful modelling. 



through the currently observed orbital period, spectral type, 
mass ratio and even individual masses A4i and M2. 

So far we have not used the additional constraints from 
other observational properties of AE Aqr. These distinguish 
between different scenarios. 



5 PROPERTIES OF CVS DESCENDING FROM 
SUPERSOFT BINARIES 

5.1 The WD spin 

We first consider the current accretion state in AE Aqr. 
The rapidly rotating WD (Pspin = 33 s) is spinning down 
on a timescale of ~ lO'^ yr, consistent with powering the 
mass expulsion in the propeller phase (Wynn et al., 1997). 
As the timescale for the spin to reach equilibrium with the 
mass transfer rate is < 10^ yr, the high mass transfer rate 
which spun up the WD must have ended only about 10^ yr 
in the past. Figure 8 compares the mass transfer rate evolu- 
tion for the four full binary models (B1-B3, B6). Clearly all 
cases predict the right M2, but all except one (B6) take too 
much time to reach that position. The spindown properties 
therefore clearly favour a large initial mass ratio for the pro- 
genitor of AE Aqr, and some fractional mass accumulation 
(here an average of 0.3) on to the WD during the supersoft 
thermal-timescale mass transfer phase. 



5.2 Abundances 

Another potential way to distinguish between the models 
lies in their different initial masses and central hydrogen 
fractions. If we combine the observations of the extreme UV 
line ratios of Civ and Nv with the apparent oversize of the 
donor for its spectral type (long orbital period), the natu- 
ral explanation is chemical evolution of the stellar interior. 
Given the complex convective history and the different ini- 
tial core masses and compositions of our four cases, we can 
construct e.g. an evolutionary sequence of characteristic sur- 
face elements as shown in Fig. 9. Despite the loss of more 
than a solar mass from the donor, solar composition is still 
maintained for a long time in sequence B2. Later however 
the ratios of ^^C to ^^C and C to N decrease by more than 
one or two orders of magnitude respectively. Vertical lines 
mark the point in the evolution where the orbital period and 
mass ratio of AE Aqr are matched: clearly this sequence pro- 
vides more than enough C depletion to explain the UV ob- 
servations of Jameson et al. (1980). Very similar results are 
found in the other cases. Bearing in mind that B2 was not 
a particularly strong case of thermal-timescale mass trans- 
fer, it seems very likely that chemical anomalies like these 
are quite common in post-thermal-timescale mass transfer 
systems. 

In order to understand the origin of this specific form 
of change in the C/N ratio, we take a closer look at the in- 
ternal structure of sequence B3 in Fig. 10. Mass transfer in 
this sequence started with a 1.8 M© donor having only 19% 
hydrogen left in the core. The left panel of Fig. 10 shows 
the internal profile of C/N shortly after filling the Roche 
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Figure 8. Time evolution of the mass transfer rate of the 
model sequences Bl (dotted), B2 (full), B3 (dash-dotted), and B6 
(dashed) shown in the previous figures. After a brief ( < 10''yr) 
phase of high mass transfer rates, the system settles towards its 
current rate. Each track reaches the current mass ratio q = 0.64 
at the position marked by an asterisk within the shaded region 
which covers the estimated range of the current mass transfer 
rate. The time axis is offset such that all tracks pass through the 
beginning of the spindown phase (marked by the spiked cross and 
assumed to start at 10~® Mq yr~^) at time t = 0. Larger values 
of the inital mass ratio qi accelerate the decline and thus shorten 
the delay between spin— up and the currently observed propeller 
state, leaving model B6 as the most suitable progenitor system. 



lobe, with mass M2 — 1.69 Mq. The initial solar composi- 
tion has been modified severly by nuclear burning and sub- 
sequent mixing in the central < IMq . Although substantial 
energy generation by CNO burning is limited to the very 
central regions (especially the convective core plateau in the 
inner ~ 0.1 Mq), the relatively long duration of the previous 
MS lifetime has allowed the nuclear reactions to change the 
composition much further out. The temperature-dependent 
CNO equilibrium ratios have been reached out to O.6M0, 
whereas the slower nuclear reactions beyond that have only 
started to reduce the amount of C. Once thermal-timescale 
mass transfer sets in however, this chemical profile remains 
frozen (everything happens within < lO'^yr, cf. Fig. 8). Ma- 
terial is removed layer by layer from the top, until finally 
an outer convection zone starts to penetrate the chemi- 
cally modified region, or the region itself is laid bare al- 
ready. The right panel of Fig. 10 illustrates this, when only 
M2 = 0.69 Mq is left. The C/N ratio in the smaU surface 
convection zone is more than a factor of 100 below solar, i.e. 
even the weakest thermal-timescale mass transfer in case 
B3 has no trouble explaining the anomalous composition of 
AE Aqr. 



Figure 9. Change of surface chemical composition in sequence 
B2: The system evolves from right to left. After an initial phase 
of constant (solar) composition, the ratios of both C/N (lower 
curve) and C-12/C-13 (upper curve, dashed) drop severely as the 
outer convective envelope reaches down into previously CNO- 
burning regions. Vertical lines mark the position where the se- 
quence reaches the orbital period (dashed) and mass ratio (dot- 
ted) of AE Aqr respectively. 
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Figure 10. Internal C/N ratio profile taken from two models 
of sequence B3. The lefthand panel shows the situation at M2 = 
1.69 Mq during thermal-timescale mass transfer, while the model 
in the righthand panel at M2 = 0.69 Mq is a close representation 
of AE Aqr now. Note the reduction of the surface ratio by a factor 
of ~ 120 once mass loss has removed the solar composition layers 
completely. 



5.3 Evolutionary endpoints: the AM CVn systems 

The later stages and endpoints of the evolutionary scenar- 
ios found for AE Aqr are connected to the advanced initial 
stage of nuclear burning on the MS. As the insets (a) in 
Figs. 6 and 7 show, the most extreme case B6 (and quite 
similar B3, which has been stopped around Porb ~ 5hr) 
reach an orbital period of 60 min or below. In this state the 
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donors consist largely of helium produced during detached 
MS burning. Thus these systems qualify as AM CVn systems 
or at least intermediate cases such as V485 Cen (Augusteijn 
et al., 1996). An almost identical channel of forming ultra- 
compact binaries has recently been studied by Podsiadlowski 
et al. (2001), focussing on initial donor masses below 1.4 M© 
and binary population synthesis aspects. 

5.4 Other properties 

(i) The upper limit for the donor mass that potentially 
leads to similar evolutionary tracks is determined by the 
occurence of a DDI, and the requirement of a strong brak- 
ing mechanism at the end of the thermal-timescale mass 
transfer supersoft phase. Both are extremely sensitive to 
the amount of mass accreted during the high M2 episode, 
which certainly needs further investigation. We suggest that 
masses up to 2Mq should be considered. 

(ii) After passing the state exemplified by AE Aqr, all 
donor stars in our models have spectral types too late for 
their orbital periods. Hence post-thermal-timescale mass 
transfer systems provide a natural explanation for the ap- 
parently evolved state of many CVs with periods above 6 hr 
(Baraffe & Kolb, 2000). Podsiadlowski et al. (2001) claim 
agreement between the properties of their population and 
the observed SpT distribution in CVs (Beuermann et al., 
1998). 

(iii) A simple argument (King & Schenker, 2002) implies 
a substantial number of evolved systems among CVs: at the 
borderline between systems brought into contact by some 
form of braking and those by nuclear evolution there will 
always be donors which have almost finished their MS phase, 
but are still kept at short orbital periods by a thermal- 
timescale mass transfer phase followed by normal angular 
momentum loss driven evolution. 

(iv) Most importantly, AE Aqr seems to be living proof 
that this kind of evolution does exist. If (as we strongly 
suspect; see below) a large number of systems with similar 
properties have passed through such a phase, the same group 
of systems may account for short-period supersoft binaries, 
AM CVns, V485 Cen-typc close binaries, and a substantial 
fraction of CVs. Although these probably include unusual 
ones like AE Aqr or V1309 Ori, many of them may appear 
fairly normal. 

(v) If we assume the track of B6 as shown in Fig. 8 as 
typical for all supersoft, AM CVn, and AE Aqr-like stages, 
one would expect roughly equal durations of the supersoft 
and spindown phases. The subsequent lifetime as normal 
CV above the period gap is about 100 times longer. More 
difficult to estimate is the fraction of post-TTMT systems 
that end up as AM CVn stars: the uncertainties of post- 
common envelope distributions, DDI, and magnetic braking 
issues prevent serious quantitative predictions. Direct com- 
parison with observed statistics is further complicated by 
poorly understood selection effects in these different classes. 

A related question concerns the expected fraction of post 
thermal-timescale mass transfer systems among CVs. We 
make a simple estimate using recently improved fits to the 
Galactic mass function (Chabricr, 2001), which enter as the 
distribution of initial donor star masses between 0.1 Mq and 
2.1 Mq. The upper limit is set by the onset of delayed dy- 
namical instability {qi = 3) for a white dwarf mass fixed for 



simplicity at 0.7 Mq. Assuming that these systems evolve 
as normal CVs for qi < 1 and via thermal-timescale mass 
transfer for qi > 1 we find that in a steady state about one- 
third of current CVs should be post-thermal-timescale mass 
transfer systems. As a simple correction for the reduced bi- 
nary fraction at low masses we have applied a weighting 
factors of 0.7 (for the thermal-timescale mass transfer mass 
range) and 0.3 (for normal CV donor masses). While this 
is obviously a fairly crude estimate, it does strongly suggest 
that a substantial fraction of CVs could show signs of an 
earlier thermal-timescale mass transfer phase. We note that 
more exhaustive studies of CV formation lead to roughly 
similar results. For example Table 1 from de Kool (1992) 
indicates that a large fraction (if not a majority) of zero- 
age binaries will undergo thermally unstable mass transfer 
rather than becoming a CV directly. 



6 CONCLUSION 

We have shown that the assumption of thermal-timescale 
mass transfer in the recent past of AE Aqr provides plausi- 
ble explanations for all of its current observational proper- 
ties. The system parameters for the best-fitting progenitor 
model presented in this paper (B6) are M2i — 1.6 Mq (fairly 
far evolved on the MS) and Pi = 18.6 hr with a primary of 
Mii = 0.6 M0 which manages to accrete about 30% of the 
total transferred mass upon reaching the current phase of 
AE Aqr. This indicates a relatively large initial mass ratio, 
and some fraction of mass accretion during the thermal- 
timescale mass transfer phase. Although the precise values 
may change, the idea is very robust, as the differential re- 
sults from Sect. 3 & 4 clearly indicate that it will be always 
possibly to end up with an excellent model for AE Aqr. 

The short duration of the thermal timescalc mass 
transfer and the transition stage (where AE Aqr currently 
is) implies a large birthrate, and thus suggests a large num- 
ber of systems passing through similar evolutions. We ex- 
pect descendants from systems similar to AE Aqr, and thus 
from supersoft binaries, to form a substantial fraction of 
the currently known CVs. The contortions of the — M2 — P 
curves of Fig. 5 suggest that non-magnetic descendant sys- 
tems may change between recurrent novae, novalike and 
dwarf nova behaviour during and after the transition from 
thermal-timescale mass transfer to normal CV evolution. 
Descendants with significant WD magnetic fields will also 
appear in various guises during these phases. We suggest 
that the long-period AM Her system V1309 Ori (Schmidt 
& Stockman, 2001; Staude et al., 2001) is a descendant of 
a supersoft binary, cf. King et al. (2002). This system is 
apparently able to synchronise at its unusually long pe- 
riod of 8 hr because of the drop in mass transfer rate at 
the end of the thermal-timescale mass transfer phase. The 
pulsing supersoft source XMMU J004319.4-f 411758 found 
by XMM in M31 may be an example of a progenitor still 
in the thermal-timescale mass transfer phase (King et al., 
2002). Note that V1309 Ori, the slightly nonsynchronous 
polar BY Cam (Bonnet-Bidaud & Mouchct, 1987) and the 
intermediate polar TX Col (Mouchet et al., 1990) all show 
high Nv/Civ ratios similar to AE Aqr. Abundance anomar 
lies possibly related to stripping of a partially evolved com- 
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panion were suggested by Mouchet et al. (1990) in the latter 
two cases. 

We conclude that in spite of its apparent uniqueness, 
AE Aqr is only the first confirmed member of a much larger 
population of post-supcrsoft binaries, constituting a signif- 
icant fraction of all CVs. A way of checking for this popula- 
tion is to determine the C/N ratio by measuring CiV 1550 
versus Nv 1238. A low value here will be strongly suggestive 
of descent from a supersoft binary. 
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